Application of Phylodynamic Tools to Inform the Public Health Response to COVID-19: Qualitative Analysis of Expert Opinions.

BACKGROUND: In the wake of the SARS-CoV-2 pandemic, scientists have scrambled to collect and analyze SARS-CoV-2 genomic data to inform public health responses to COVID-19 in real time. Open source phylogenetic and data visualization platforms for monitoring SARS-CoV-2 genomic epidemiology have rapidly gained popularity for their ability to illuminate spatial-temporal transmission patterns worldwide. However, the utility of such tools to inform public health decision-making for COVID-19 in real time remains to be explored. OBJECTIVE: The aim of this study is to convene experts in public health, infectious diseases, virology, and bioinformatics-many of whom were actively engaged in the COVID-19 response-to discuss and report on the application of phylodynamic tools to inform pandemic responses. METHODS: In total, 4 focus groups (FGs) occurred between June 2020 and June 2021, covering both the pre- and postvariant strain emergence and vaccination eras of the ongoing COVID-19 crisis. Participants included national and international academic and government researchers, clinicians, public health practitioners, and other stakeholders recruited through purposive and convenience sampling by the study team. Open-ended questions were developed to prompt discussion. FGs I and II concentrated on phylodynamics for the public health practitioner, while FGs III and IV discussed the methodological nuances of phylodynamic inference. Two FGs per topic area to increase data saturation. An iterative, thematic qualitative framework was used for data analysis. RESULTS: We invited 41 experts to the FGs, and 23 (56%) agreed to participate. Across all the FG sessions, 15 (65%) of the participants were female, 17 (74%) were White, and 5 (22%) were Black. Participants were described as molecular epidemiologists (MEs; n=9, 39%), clinician-researchers (n=3, 13%), infectious disease experts (IDs; n=4, 17%), and public health professionals at the local (PHs; n=4, 17%), state (n=2, 9%), and federal (n=1, 4%) levels. They represented multiple countries in Europe, the United States, and the Caribbean. Nine major themes arose from the discussions: (1) translational/implementation science, (2) precision public health, (3) fundamental unknowns, (4) proper scientific communication, (5) methods of epidemiological investigation, (6) sampling bias, (7) interoperability standards, (8) academic/public health partnerships, and (9) resources. Collectively, participants felt that successful uptake of phylodynamic tools to inform the public health response relies on the strength of academic and public health partnerships. They called for interoperability standards in sequence data sharing, urged careful reporting to prevent misinterpretations, imagined that public health responses could be tailored to specific variants, and cited resource issues that would need to be addressed by policy makers in future outbreaks. CONCLUSIONS: This study is the first to detail the viewpoints of public health practitioners and molecular epidemiology experts on the use of viral genomic data to inform the response to the COVID-19 pandemic. The data gathered during this study provide important information from experts to help streamline the functionality and use of phylodynamic tools for pandemic responses.

Persistence of SARS-CoV-2 omicron variant in children and utility of rapid antigen testing as an indicator of culturable virus.

We screened 65 longitudinally-collected nasal swab samples from 31 children aged 0-16 years who were positive for SARS-CoV-2 omicron BA.1. By day 7 after onset of symptoms 48% of children remained positive by rapid antigen test. In a sample subset we found 100% correlation between antigen test results and virus culture.

Structural Evolution of Delta (B.1.617.2) and Omicron (BA.1) Spike Glycoproteins.

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627-635 region, on a small region of the receptor binding motif (aa 483-485), and the so-called "fusion-peptide proximal region". In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.

Severe Acute Respiratory Syndrome Coronavirus 2 Delta Vaccine Breakthrough Transmissibility in Alachua County, Florida.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant has caused a dramatic resurgence in infections in the United Sates, raising questions regarding potential transmissibility among vaccinated individuals. METHODS: Between October 2020 and July 2021, we sequenced 4439 SARS-CoV-2 full genomes, 23% of all known infections in Alachua County, Florida, including 109 vaccine breakthrough cases. Univariate and multivariate regression analyses were conducted to evaluate associations between viral RNA burden and patient characteristics. Contact tracing and phylogenetic analysis were used to investigate direct transmissions involving vaccinated individuals. RESULTS: The majority of breakthrough sequences with lineage assignment were classified as Delta variants (74.6%) and occurred, on average, about 3 months (104 ±â€…57.5 days) after full vaccination, at the same time (June-July 2021) of Delta variant exponential spread within the county. Six Delta variant transmission pairs between fully vaccinated individuals were identified through contact tracing, 3 of which were confirmed by phylogenetic analysis. Delta breakthroughs exhibited broad viral RNA copy number values during acute infection (interquartile range, 1.2-8.64 Log copies/mL), on average 38% lower than matched unvaccinated patients (3.29-10.81 Log copies/mL, P < .00001). Nevertheless, 49% to 50% of all breakthroughs, and 56% to 60% of Delta-infected breakthroughs exhibited viral RNA levels above the transmissibility threshold (4 Log copies/mL) irrespective of time after vaccination. CONCLUSIONS: Delta infection transmissibility and general viral RNA quantification patterns in vaccinated individuals suggest limited levels of sterilizing immunity that need to be considered by public health policies. In particular, ongoing evaluation of vaccine boosters should specifically address whether extra vaccine doses curb breakthrough contribution to epidemic spread.

Low-frequency variants in mildly symptomatic vaccine breakthrough infections presents a doubled-edged sword.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) has raised questions regarding vaccine protection against SARS-CoV-2 infection, transmission, and ongoing virus evolution. Twenty-three mildly symptomatic "vaccination breakthrough" infections were identified as early as January 2021 in Alachua County, Florida, among individuals fully vaccinated with either the BNT162b2 (Pfizer) or the Ad26 (Janssen/J&J) vaccines. SARS-CoV-2 genomes were successfully generated for 11 of the vaccine breakthroughs, and 878 individuals in the surrounding area and were included for reference-based phylogenetic investigation. These 11 individuals were characterized by infection with VOCs, but also low-frequency variants present within the surrounding population. Low-frequency mutations were observed, which have been more recently identified as mutations of interest owing to their location within targeted immune epitopes (P812L) and association with increased replicative capacity (L18F). We present these results to posit the nature of the efficacy of vaccines in reducing symptoms as both a blessing and a curse-as vaccination becomes more widespread and self-motivated testing reduced owing to the absence of severe symptoms, we face the challenge of early recognition of novel mutations of potential concern. This case study highlights the critical need for continued testing and monitoring of infection and transmission among individuals regardless of vaccination status.

Rapid Emergence and Spread of Severe Acute Respiratory Syndrome Coronavirus 2 Gamma (P.1) Variant in Haiti.

After an initial wave of coronavirus disease 2019 (COVID-19) in Haiti in summer 2020 (primarily lineage B.1), seropositivity for anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) immunoglobulin G (IgG) was ~40%. Variant P.1 (gamma) was introduced in February 2021, with an initially limited introduction followed by exponential local dissemination within this unvaccinated population with prior exposure to earlier SARS-CoV-2 lineages.

Rapid SARS-CoV-2 diagnosis using disposable strips and a metal-oxide-semiconductor field-effect transistor platform.

The SARS-CoV-2 pandemic has had a significant impact worldwide. Currently, the most common detection methods for the virus are polymerase chain reaction (PCR) and lateral flow tests. PCR takes more than an hour to obtain the results and lateral flow tests have difficulty with detecting the virus at low concentrations. In this study, 60 clinical human saliva samples, which included 30 positive and 30 negative samples confirmed with RT-PCR, were screened for COVID-19 using disposable glucose biosensor strips and a reusable printed circuit board. The disposable strips were gold plated and functionalized to immobilize antibodies on the gold film. After functionalization, the strips were connected to the gate electrode of a metal-oxide-semiconductor field-effect transistor on the printed circuit board to amplify the test signals. A synchronous double-pulsed bias voltage was applied to the drain of the transistor and strips. The resulting change in drain waveforms was converted to digital readings. The RT-PCR-confirmed saliva samples were tested again using quantitative PCR (RT-qPCR) to determine cycling threshold (Ct) values. Ct values up to 45 refer to the number of amplification cycles needed to detect the presence of the virus. These PCR results were compared with digital readings from the sensor to better evaluate the sensor technology. The results indicate that the samples with a range of Ct values from 17.8 to 35 can be differentiated, which highlights the increased sensitivity of this sensor technology. This research exhibits the potential of this biosensor technology to be further developed into a cost-effective, point-of-care, and portable rapid detection method for SARS-CoV-2.

Independent infections of porcine deltacoronavirus among Haitian children.

Coronaviruses have caused three major epidemics since 2003, including the ongoing SARS-CoV-2 pandemic. In each case, the emergence of coronavirus in our species has been associated with zoonotic transmissions from animal reservoirs1,2, underscoring how prone such pathogens are to spill over and adapt to new species. Among the four recognized genera of the family Coronaviridae, human infections reported so far have been limited to alphacoronaviruses and betacoronaviruses3-5. Here we identify porcine deltacoronavirus strains in plasma samples of three Haitian children with acute undifferentiated febrile illness. Genomic and evolutionary analyses reveal that human infections were the result of at least two independent zoonoses of distinct viral lineages that acquired the same mutational signature in the genes encoding Nsp15 and the spike glycoprotein. In particular, structural analysis predicts that one of the changes in the spike S1 subunit, which contains the receptor-binding domain, may affect the flexibility of the protein and its binding to the host cell receptor. Our findings highlight the potential for evolutionary change and adaptation leading to human infections by coronaviruses outside of the previously recognized human-associated coronavirus groups, particularly in settings where there may be close human-animal contact.

Isolation of a Novel Recombinant Canine Coronavirus From a Visitor to Haiti: Further Evidence of Transmission of Coronaviruses of Zoonotic Origin to Humans.

We isolated a novel coronavirus from a medical team member presenting with fever and malaise after travel to Haiti. The virus showed 99.4% similarity with a recombinant canine coronavirus recently identified in a pneumonia patient in Malaysia, suggesting that infection with this virus and/or recombinant variants occurs in multiple locations.

Comparison of SARS-CoV-2 Receptors Expression in Primary Endothelial Cells and Retinoic Acid-Differentiated Human Neuronal Cells.

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is primarily responsible for coronavirus disease (COVID-19) and it is characterized by respiratory illness with fever and dyspnea. Severe vascular problems and several other manifestations, including neurological ones, have also been frequently reported, particularly in the great majority of "long hauler" patients. SARS-CoV-2 infects and replicates in lung epithelial cells, while dysfunction of endothelial and neuronal brain cells has been observed in the absence of productive infection. It has been shown that the Spike protein can interact with specific cellular receptors, supporting both viral entry and cellular dysfunction. It is thus clear that understanding how and when these receptors are regulated, as well as how much they are expressed would help in unveiling the multifaceted aspects of this disease. Here, we show that SH-SY5Y neuroblastoma cells express three important cellular surface molecules that interact with the Spike protein, namely ACE2, TMPRSS2, and NRP1. Their levels increase when cells are treated with retinoic acid (RA), a commonly used agent known to promote differentiation. This increase matched the higher levels of receptors observed on HUVEC (primary human umbilical vein endothelial cells). We also show by confocal imaging that replication-defective pseudoviruses carrying the SARS-CoV-2 Spike protein can infect differentiated and undifferentiated SH-SY5Y, and HUVEC cells, although with different efficiencies. Neuronal cells and endothelial cells are potential targets for SARS-CoV-2 infection and the interaction of the Spike viral protein with these cells may cause their dysregulation. Characterizing RNA and protein expression tempo, mode, and levels of different SARS-CoV-2 receptors on both cell subpopulations may have clinical relevance for the diagnosis and treatment of COVID-19-infected subjects, including long hauler patients with neurological manifestations.

Collection of SARS-CoV-2 Virus from the Air of a Clinic Within a University Student Health Care Center and Analyses of the Viral Genomic Sequence.

The progression of COVID-19 worldwide can be tracked by identifying mutations within the genomic sequence of SARS-CoV-2 that occur as a function of time. Such efforts currently rely on sequencing the genome of SARS-CoV-2 in patient specimens (direct sequencing) or of virus isolated from patient specimens in cell cultures. A pilot SARS-CoV-2 air sampling study conducted at a clinic within a university student health care center detected the virus vRNA, with an estimated concentration of 0.87 virus genomes L-1 air. To determine whether the virus detected was viable ('live'), attempts were made to isolate the virus in cell cultures. Virus-induced cytopathic effects (CPE) were observed within two days post-inoculation of Vero E6 cells with collection media from air samples; however, rtRT-PCR tests for SARS-CoV-2 vRNA from cell culture were negative. Instead, three other fast-growing human respiratory viruses were isolated and subsequently identified, illustrating the challenge in isolating SARS-CoV-2 when multiple viruses are present in a test sample. The complete SAR-CoV-2 genomic sequence was nevertheless determined by Sanger sequencing and most closely resembles SARS-CoV-2 genomes previously described in Georgia, USA. Results of this study illustrate the feasibility of tracking progression of the COVID-19 pandemic using environmental aerosol samples instead of human specimens. Collection of a positive sample from a distance more than 2 m away from the nearest patient traffic implies the virus was in an aerosol.

Multiple Recombination Events and Strong Purifying Selection at the Origin of SARS-CoV-2 Spike Glycoprotein Increased Correlated Dynamic Movements.

Our evolutionary and structural analyses revealed that the severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) spike gene is a complex mosaic resulting from several recombination events. Additionally, the fixation of variants has mainly been driven by purifying selection, suggesting the presence of conserved structural features. Our dynamic simulations identified two main long-range covariant dynamic movements of the novel glycoprotein, and showed that, as a result of the evolutionary duality, they are preserved. The first movement involves the receptor binding domain with the N-terminal domain and the C-terminal domain 2 and is maintained across human, bat and pangolin coronaviruses. The second is a complex network of long-range dynamics specific to SARS-CoV-2 involving the novel PRRA and the conserved KR*SF cleavage sites, as well as conserved segments in C-terminal domain 3. These movements, essential for host cell binding, are maintained by hinges conserved across human, bat, and pangolin coronaviruses glycoproteins. The hinges, located around Threonine 333 and Proline 527 within the N-terminal domain and C-terminal domain 2, represent candidate targets for the future development of novel pan-coronavirus inhibitors. In summary, we show that while recombination created a new configuration that increased the covariant dynamic movements of the SARS-CoV-2 glycoprotein, negative selection preserved its inter-domain structure throughout evolution in different hosts and inter-species transmissions.

Correction: A Snapshot of SARS-CoV-2 Genome Availability up to April 2020 and its Implications: Data Analysis.

[This corrects the article DOI: 10.2196/19170.].

Multiple Recombination Events and Strong Purifying Selection at the Origin of SARS-CoV-2 Spike Glycoprotein Increased Correlated Dynamic Movements

Our evolutionary and structural analyses revealed that the severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) spike gene is a complex mosaic resulting from several recombination events. Additionally, the fixation of variants has mainly been driven by purifying selection, suggesting the presence of conserved structural features. Our dynamic simulations identified two main long-range covariant dynamic movements of the novel glycoprotein, and showed that, as a result of the evolutionary duality, they are preserved. The first movement involves the receptor binding domain with the N-terminal domain and the C-terminal domain 2 and is maintained across human, bat and pangolin coronaviruses. The second is a complex network of long-range dynamics specific to SARS-CoV-2 involving the novel PRRA and the conserved KR*SF cleavage sites, as well as conserved segments in C-terminal domain 3. These movements, essential for host cell binding, are maintained by hinges conserved across human, bat, and pangolin coronaviruses glycoproteins. The hinges, located around Threonine 333 and Proline 527 within the N-terminal domain and C-terminal domain 2, represent candidate targets for the future development of novel pan-coronavirus inhibitors. In summary, we show that while recombination created a new configuration that increased the covariant dynamic movements of the SARS-CoV-2 glycoprotein, negative selection preserved its inter-domain structure throughout evolution in different hosts and inter-species transmissions.

Collection of SARS-CoV-2 Virus from the Air of a Clinic within a University Student Health Care Center and Analyses of the Viral Genomic Sequence

The progression of COVID-19 worldwide can be tracked by identifying mutations within the genomic sequence of SARS-CoV-2 that occur as a function of time. Such efforts currently rely on sequencing the genome of SARS-CoV-2 in patient specimens (direct sequencing) or of virus isolated from patient specimens in cell cultures. A pilot SARS-CoV-2 air sampling study conducted at a clinic within a university student health care center detected the virus vRNA, with an estimated concentration of 0.87 virus genomes L-1 air. To determine whether the virus detected was viable ('live'), attempts were made to isolate the virus in cell cultures. Virus-induced cytopathic effects (CPE) were observed within two days post-inoculation of Vero E6 cells with collection media from air samples;however, rtRT-PCR tests for SARS-CoV-2 vRNA from cell culture were negative. Instead, three other fast-growing human respiratory viruses were isolated and subsequently identified. illustrating the challenge in isolating SARS-CoV-2 when multiple viruses are present in a test sample. The complete SAR-CoV-2 genomic sequence was nevertheless determined by Sanger sequencing and most closely resembles SARS-CoV-2 genomes previously described in Georgia, USA. Results of this study illustrate the feasibility of tracking progression of the COVID-19 pandemic using environmental aerosol samples instead of human specimens. Collection of a positive sample from a distance more than 2 m away from the nearest patient traffic implies the virus was in an aerosol.

A snapshot of SARS-CoV-2 genome availability up to April 2020 and its implications

BACKGROUND: : The SARS-CoV-2 pandemic has been growing exponentially, affecting over four million people and causing enormous distress to economies and societies worldwide. A plethora of analyses based on viral sequences has already been published in scientific journals, as well as through non-peer reviewed channels, to investigate SARS-CoV-2 genetic heterogeneity and spatiotemporal dissemination. Yet, a systematic investigation of phylogenetic information and sampling bias in the available data is missing. OBJECTIVE: The objective of this study was to determine the quality of the current SARS-CoV-2 full genome data, in terms of sampling bias as well as phylogenetic and temporal signal, to inform and guide the scientific community. METHODS: We used maximum likelihood based methods to assess the presence of sufficient information for robust phylogenetic and phylogeographic studies in several SARS-CoV-2 sequence alignments, assembled from GISAID data released between March and April 2020. RESULTS: Although number of high quality full-genomes is growing daily, and recent sequence data released in April contains sufficient phylogenetic information that would allow reliable inference of phylogenetic relationships, country-specific SARS-CoV-2 data sets still present severe limitations. CONCLUSIONS: At the present time, studies assessing within country spread or transmission clusters should be considered preliminary at best, or hypothesis generating. Hence, the need for interpreting current reports with caution, and continuing concerted efforts to increase number and quality of sequences required for robust tracing of the epidemic. CLINICALTRIAL: INTERNATIONAL REGISTERED REPORT: RR2-https://doi.org/10.1101/2020.04.01.020594.